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Abstract 

We investigate the statistics of eigenstates in a weak self-affine disordered po- 
tential in one dimension, whose Gaussian fluctuations grow with distance with a 
positive Hurst exponent H. Typical eigenstates are superlocalized on samples much 
larger than a well-defined crossover length, which diverges in the weak-disorder 
regime. We present a parallel analytical investigation of the statistics of these su- 
perlocalized states in the discrete and the continuum formalisms. For the discrete 
tight-binding model, the effective localization length decays logarithmically with 
the sample size, and the logarithm of the transmission is marginally self-averaging. 
For the continuum Schrodinger equation, the superlocalization phenomenon has 
more drastic effects. The effective localization length decays as a power of the 
sample length, and the logarithm of the transmission is fully non-self-averaging. 
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1 Introduction 



The Anderson localization of a quantum-mechanical particle by a random potential is 
by now an old and well-understood problem [1]. This is especially so in one dimension. 
Consider for definiteness the tight-binding equation: 

Ipn+l + 1pn-l + Vnlpn ^ Elpn- (1-1) 

In the usual situation where the site potentials are uncorrelated, all the eigenstates 
are known to be exponentially locahzed (see e.g. [2, 3]). 

The peculiar features of the localization problem for various kinds of one- dimensional 
disordered potentials with non-trivial correlations have also been investigated, includ- 
ing dimer models [4], potentials with power-law correlations [5], potentials generated by 
chaotic dynamical systems [6], and potentials whose correlations are designed on pur- 
pose [7] . These examples share the common feature of spatial stationarity: the statistics 
of the potential is invariant under translation, so that its two-point correlation has a 
well-defined thermodynamical limit (KnV^) = C^_„, which only depends on the distance 
\m — n\. 

The more exotic situation of non- stationary random potentials, whose fluctuations 
grow with distance, has only been considered more recently [8, 9, 10, 11, 12, 13]. An 
example of most physical interest is that of a self-affine Gaussian potential with Hurst 
exponent < if < 1, such that 

((Kn-K)')- A2|m-n|2^. (1.2) 

Such a sequence of potentials can be generated by fractional Brownian motion, with a 
proper choice of stationary but correlated increments 

en = Vn- Vn-1. (1.3) 

The particular case of the usual Brownian motion, corresponding to stationary and in- 
dependent increments ((ern^n) — ^^mn)-i has a Hurst exponent H = 1/2. 

An alternative way of characterizing non-stationary potentials consists in considering 
their structure factor S{q) = {V{q)V{—q)) , assuming a power-law divergence of the form 

S(q) - \qr (1-4) 

in the long- wavelength limit (g ^ 0). Only the smaller values of the scaling exponent 
(a < 1) correspond to stationary potentials, with long-range correlations falling off as 
Cn ~ \n\^^^~°'\ Larger values of a necessarily correspond to non-stationary potentials. 
The above self-affine potentials with stationary increments are obtained in the range 
1 < a < 3, with the correspondence a = 2H + 1. 

The main novel feature induced by the non- stationarity of the potential is that the 
effective disorder strength depends on the spatial scale. The typical potential fluctuation 
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over a distance indeed grows as Kms(-^) — AA^^. It becomes therefore comparable 
with the bandwidth when the distance reaches the crossover length 

TV,, = A-V^. (1.5) 

The situation of most interest is that of a weak disorder (A ^ 1), so that A^^cr ^ 1- 
The eigenstates are conventionally localized on scales N <^ N^, with a localization 
length scaling as ^ ~ 1/A^. On larger scales {N » A^cr), however, eigenstates become 
strongly localized or superlocalized [12, 13], because wavefunctions fall off very fast in 
the classically forbidden zones of the potential {\E — Vn\ > 2). Let us mention that 
an alternative viewpoint [8, 9] consists in keeping fixed the effective potential width 
K'ms(^) = AA^^ = S at the scale of the system size A^. The microscopic disorder 
strength A = H/N^ therefore gets rescaled to smaller and smaller values. The crossover 
length A^cr = N/T}/^ grows proportionally to the sample size, whereas the localization 
length scales as ^ ~ N"^^ . For a large enough non-stationarity {H > 1/2), such that ^ may 
become much larger than the sample size N, a crossover line in the (E, E) plane between 
extended and localized states has been observed [8, 9], and theoretically explained [12. 13]. 

The goal of the present work is to provide the first quantitative analysis of the statis- 
tics of superlocalized eigenstates, especially in the regime A^ ^ A'cr, where the super- 
localization phenomenon is fully developed. We shall successively deal with the discrete 
tight- binding model (Section 2) and the continuum Schrodinger equation (Section 3). 

The key quantity considered throughout this work is the effective Lyapunov exponent 

7iv = -^InlV'jvl, (1.6) 

where ^j)n is the generic (growing) solution to (1.1). The effective Lyapunov exponent 
provides an estimate of the global growth rate of this solution over N lattice sites, and 
therefore of the effective decay rate of eigenstates over the same range. In other words, 
the effective localization length at scale N is 

iN = — . (1.7) 

In 

The effective Lyapunov exponent is also a central quantity in the theory of coherent 
quantum transport. Indeed the two-probe Landauer formula [14] expresses the zero- 
temperature conductance oia. sample made of A^ lattice sites, in terms of the intensity 
transmission T/v across that sample, as 

2e^ 

Qn^-^Tm. (1.8) 

Furthermore, in the insulating regime where the transmission is small, it is known to 
scale as T/v ~ I/I'V^at^, hence 

InT^v ~ -2N'^N- (1.9) 

The effective Lyapunov exponent, and related physical quantities such as the effective 
localization length and the logarithm of the transmission (conductance), exhibit different 
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kinds of scaling behavior in the two situations to be investigated successively, the discrete 
tight-binding model (Section 2) and the continuum Schrodinger equation (Section 3). 
This basic difference is further commented on in the Discussion (Section 4). 



2 The discrete tight-binding model 
2.1 Reminder 

We find it useful to first give a brief reminder on the conventional situation of a stationary 
random potential. The statistics of the size-dependent effective Lyapunov exponent 7Ar is 
then very similar to that of the free-energy density of a disordered statistical-mechanical 
system [3, 15, 16]. The effective Lyapunov exponent has a well-defined self-averaging 
limit 7 in the N ^ oo limit, simply referred to as the Lyapunov exponent. Its reciprocal 
^ = 1/7 is interpreted as the localization length of the problem. 

The fluctuations of 7^^ around 7 are Gaussian and scale as var7jv = {{'Jn — lY) ~ 
a"^ /N for a finite but large enough sample (A'"7 ^ 1). More generally, the product A'"7;v 
is extensive, in the strong sense that all its cumulants 

{{{NinT)) ~ CfeTV (2.1) 

grow linearly with N [3], with amplitudes Ci = 7, C2 = cr^, and so on. It is worth noticing 
that N'~fi\f, which is analogous to the total free energy, is also a physical observable in the 
present context, as it is nothing but the logarithm of the transmission (conductance), up 
to a constant factor [see (1.9)]. 

Furthermore, in the usual situation of independent site potentials with (V^) = and 
{VmVn) = W^^mn, fhc Lyapuuov exponent behaves as follows in the weak-disorder limit 
{W^ < 1): 



Inside the band, i.e., for l^l < 2, setting E = 2cosg, the celebrated result [2] 



Ssin^g ~ 2(4 -£;^) 



7 ^ ^7— 2- = (2-2) 



shows that the locahzation length diverges as l/W"^. 
Outside the band, i.e., for > 2, setting \E\ = 2 cosh t, i.e., 

.^ln l^l + f^ >0. (2.3) 

the generic solution to (1.1) grows exponentially as ~ e"* in the absence of 
disordered potential, hence 

7^t iW'^^Q). (2.4) 
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• Near band edges, i.e., for £■ — > ±2 and W'^ — > simultaneously, the Lyapunov 
exponent obeys a scaling law of the form 

where the scaling function F is explicitly known in terms of Airy functions [17]. 
2.2 General results 

We now turn to the current problem, namely the tight-binding equation (1.1) with a 
weak self-afRne potential whose fluctuations grow with distance according to (1.2), with 

<C 1. We choose the zero of energies as the site potential at the origin (Vq = 0), and 
assume that the origin sits in an allowed zone, i.e., \E\ < 2. 

A long enough sample {N » A^cr) generically consists of an alternation of allowed 
zones {\E — Vn\ < 2) and of forbidden zones {\E — Vn\ > 2). It can be argued from 
the above results that the growth rate of a generic wavefunction is proportional to 
in allowed zones, i.e., inside the 'local band' around site n (where energies are shifted 
by Vn) [see (2.2)], whereas it is of order unity in forbidden zones, i.e., outside the 'local 
band' [see (2.4)]. 

This picture is fully corroborated by the plots shown in References [12, 13], where 
eigenstates are indeed seen to be essentially constant in allowed zones, and to drop 
very suddenly in forbidden zones. Our main goal is to turn the above picture into a 
quantitative description of the statistics of the effective Lyapunov exponent (1.6). 

Our starting point is the approximate formula 

ipN ^ exp (^tn\ , (2.6) 

\n=0 / 



with the definition 
hence 



E -Vn^2coshtn, (2.7) 



N-l 



7iv = ^ E R-etn, (2.8) 

n=0 




^ > in forbidden zones {\E-Vn\>2), 



with 



Re tr., 

in allowed zones {\E — Vn\ < 2). 

(2.9) 

Let us first comment on (2.6) for a while. This expression is a full discrete analogue of 
the W.K.B. integral (3.4), to be used in Section 3. It provides a quantitative description 
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of the growing solution to (1.1) in the forbidden zones of the potential, under the sole 
hypothesis that the sequence of site potentials has small increments: 



en = Vn- K-i < 1. (2.10) 
Expression (2.6) can be easily derived by introducing the Riccati variables [16, 18, 19] 

Rn^^. (2.11) 

Assuming ■^q = 1; we have ipn = Rq ■ ■ ■ Rn-i, and 

N-l 

N 



7iv = 4Eln|^n|. (2.12) 



n=0 

The variables Rn obey the recursion 



If we now set 



the variables Xn obey 



with 



Rn^E-Vn- -J—. (2.13) 



'^'^ ~ 1 r> ' R'^ ~ 1 t ' (2-14) 

1 — e"Rr,. 1 — e*"Xr. 



-2t„ + 



Xn = e-^*" ' , (2.15) 

1 + OnXn-l 



""-e*n+tn-i_i -4sinh2i^' (2.16) 

to leading order as £„ <^ 1. To the same order, the recursion (2.15) can be hnearized 
to Xn ~ e~^*"(x„_i + 5n), implying that the x„'s are proportional to the £„'s. One 
has therefore i?„ exp(t„), up to terms proportional to the £„'s. The expressions (2.6) 
and (2.8) follow at once. Keeping track of higher powers of the e„'s in the above equations 
would be an efficient way of performing systematic weak-disorder expansions [16, 19]. 

Our analysis of the statistics of the effective Lyapunov exponent 7jv is based on the 
estimate (2.8). It turns out to be advantageous to introduce the Laplace representation 

/ds 
- — cosh(s£;)Xo(2s) exp(sK), (2.17) 
ITTS 

where Kq is the modified Bessel function, and the integration contour is a vertical line 
with Res > 0. Using the Gaussian statistics of the potentials V^, with Vq = and (1.2), 
we have 

(exp(sK)) = cxp (lA^sV^) , (2.18) 

and similar expressions for higher-order characteristic functions such as {cxp{siVm + 
S2Vn))- We can therefore express the correlation functions of the variables Ret^ as 
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multiple contour integrals. By means of (2.8), the moments of the effective Lyapunov 
exponent can, in turn, be expressed as multiple integrals. The rescaling x = m/Ncr, 
y = n/Ncr implies that, as could be expected, the final results only depend on the system 
size through the dimensionless ratio 

X = — = A^/^AT, (2.19) 

where the crossover scale A^cr has been introduced in (1.5). We thus obtain 

(7^) [""dx cosh{sE) Ko{2s) exp ( ^s'^x''") , (2.20) 

X Jo J ITTS ^ ' 

(tIt) r dx dy cosh(siE) Ko{2s,) [ ^ cosh(s2^) i^o(2s2) 

X"^ Jo Jo J ITTSi J mS2 

X exp + s^s,ix''' + y^^ - \x - y^) + sly'"")) , (2.21) 

and so on. 

The regimes of short samples (A^ <^ A^cr) and of long samples (A^ 3> A^cr) deserve to 
be considered separately. 



2.3 Short samples 

We first consider relatively short samples, such that N -C A'cr, i-e., X <^ 1. In this 
regime, and for |£'| < 2, the contour integral in (2.20) is dominated by a saddle point at 
Sc ~ (2 — \E\)/x'^^ ^ 1. We thus obtain the following exponentially small estimate 

(7Ar)-exp p^JfJ^' ) (2.22) 

for the mean effective Lyapunov exponent. 

The right-hand side of (2.22) has a simple interpretation. It scales indeed as the 
probability that E — Vn has reached the closest forbidden zone, i.e., E — Vn = 2 if E > 
and £^ — Vat = — 2 if £■ < 0, so that the superlocalization phenomenon just sets in. 

For a weak but finite disorder strength A, the mean effective Lyapunov exponent {'Jn) 
is already of order A^ in the regime of usual localization, i.e., for X -C 1. The onset 
of superlocalization manifests itself as a crossover of {'Jn) fo the behavior (2.22), and it 
takes place for a system size such that 

X~|lnA|-V(2^)«l, i.e., N^(A\lnA\'/^y'^" . (2.23) 



2.4 Long samples 

The superlocalization phenomenon is fully developed in the converse regime of long sam- 
ples (A^ S> A'cr, i.e., X 3> 1). The contour integral in (2.20) is then dominated by the 
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branch cut of Kq{2s) !v —(Ins + C) at s = 0, with C being Euler's constant. We thus 
obtain 

(7iv) ^ H{lnX - 1) - £±^. (2.24) 

The mean effective Lyapunov exponent thus diverges logarithmically with the system 
size N ior N Na, irrespective of the energy E. 

As a matter of fact, the full distribution of the effective Lyapunov exponent simplifies 
in the limit of a long sample. Indeed, for n ^ A^cr one has (V^) ^ 1, so that most often 
\Vn\ ^ 1, hence Ret„ ~ In |V^|, and 

^ N-l 

7iV~ T7 E^^l^nl- (2.25) 

n=0 

The rescaling x = n/N yields the following asymptotic form of the effective Lyapunov 
exponent of long samples: 

^ HlnX + Gh ^ H\n{N/N^,) + Gh^ \n{AN") + Gh, (2.26) 

where the random additive term Gh is the following non-linear functional of the normal- 
ized fractional Brownian motion: 

Gh^ rdxln\v{x)\. (2.27) 
Jo 

The expression (2.26), which holds irrespectively of the energy E, is the main result of this 
section. The functional Gh is investigated in the Appendix. It has a universal distribution 
which only depends on the Hurst exponent H. In the case H — 1/2, corresponding to 
usual Brownian motion, similar non-linear integral functionals have been met in several 
physical problems [20]. To our knowledge, the distribution of Gh is not known explicitly, 
even in the Brownian case. Expressions for the mean and the variance of Gh are derived in 
the Appendix for all H, as well as the full distribution of Gi in the ballistic limit {H = 1). 
The mean {Gh) has a simple hnear dependence (A. 6) on the Hurst exponent H, which 
can already be read off from (2.24), while the more complex expression (A. 17) for the 
variance yqxGh only simplifies for if = 0, 1/2, and 1. Figure 1 illustrates the dependence 
of var Gh on the exponent H, whereas Figure 2 shows a plot of the probability density 
of Gi (ballistic case: exact expression (A. 33)) and of G1/2 (Brownian case), the latter 
being measured from an ensemble of numerically generated, and suitably rescaled, long 
random walks. The distributions is observed to be clearly asymmetric (skew), with very 
fast decaying tails at large values of Gh, and at small values oi Gh for any H < 1. 

2.5 Numerical results 

We have confronted our analytical predictions with numerical results in the Brownian case 
{H — 1/2). This situation is indeed simpler to handle than the generic one. On the one 
hand, the integrals over the rescaled spatial positions x, y in the predictions (2.20), (2.21) 
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Figure 1: Plot of the variance of the functional Gh entering the result (2.26), evaluated from 
the exact expression (A. 17), against the Hurst exponent H. Symbols: results (A. 25), (A. 27), 
and (A. 30) for the particular cases H = 0, 1/2, and 1. 

are elementary. On the other hand, sequences of Brownian potentials Vn can easily be 
generated numerically, by summing independent increments = ±A with A ^ 1, 
according to (1.3), with Vq = E = 0. 

The growing solution ipn has been evaluated by means of the recursion (2.13) for the 
Riccati variables with a random initial condition {Ri = tan0 with a uniform angle 0). 
The effective Lyapunov exponent has been measured by means of (2.12). The mean and 
the variance of the effective Lyapunov exponent are respectively shown in Figures 3 
and 4. The continuous curves represent the analytical results (2.20), (2.21), whereas 
symbols show numerical data for 10^ samples with two strengths of disorder, A = 0.1 
and A = 0.2. The accuracy of the data collapse and of the quantitative agreement with 
theoretical predictions provides a convincing check of our analysis. 

3 The continuum Schrodinger equation 

We now turn to the Schrodinger equation in a one- dimensional potential V{x), which 
reads 

-0"(x) + V{x)ilj{x) = ExIj{x), (3.1) 

in reduced units. The potential V{x) is again assumed to be a weak self-affine Gaussian 
disordered potential with Hurst exponent < H < 1, such that 

mx)-V{y)r)^A'\x-y\'", (3.2) 

with V{0) = and A^ < 1. 
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Figure 2: Plot of the probability density of the functional Gr entering the result (2.26). Thin 
full line: exact result (A.33) for = 1 (ballistic limit). Line with symbols (showing histogram 
bins): numerical result for = 1/2 (Brownian case). 



3.1 General results 

The present situation shares most of the characteristic features of the tight-binding prob- 
lem investigated in Section 2. Its analysis will therefore only be described succinctly. 

A long enough sample is again generically an alternation of classically allowed zones 
{V{x) < E) and of forbidden zones {y{x) > E). We assume £■ > 0, so that the origin 
belongs to an allowed zone. The growth rate of the generic solution of (3.1) up to 
distance L is measured by the effective Lyapunov exponent 

7(L) = iln|^(L)|, (3.3) 

The local growth rate of wavefunctions is again proportional to in allowed zones, 
whereas it is of order unity in forbidden zones. It is therefore legitimate to use the 
celebrated W.K.B. integral (see e.g. [21, 22]) 

^{x) ~ exp dy ^{V{y) - £;)+) , (3.4) 



with the definition 



iy{y) - E). 



V{y) — E in forbidden zones {V{y) > E), 



. in allowed zones {Viv) < E). 



(3.5) 



The integral formula (3.4) provides a quantitative description, with exponential ac- 
curacy, of the growing solution to (3.1) in the forbidden zones of the potential, whenever 
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In(NA^) 

Figure 3: Plot of the mean effective Lyapunov exponent for a Brownian potential (H = 1/2) 
with E = 0, against the logarithm of the scaling variable X = iVA^. Full curve: analytical 
prediction (2.20). Dashed line: asymptotic behavior (2.24), i.e., {'Jn) ^ (In(iVA^) — 1 — C — 
ln2)/2. Symbols: numerical data for A = 0.1 and A = 0.2. 



the latter is slowly varying. One has therefore 

liL) ^ \ dx ^{V{x) - EU. (3.6) 

Equations (3.6) and (3.5) are the continuum analogues of (2.8) and (2.9). They 
constitute the starting point of the subsequent analysis of the statistics of the effective 
Lyapunov exponent. It is again advantageous to introduce a Laplace representation: 

^{V{x)-EU = J^\[^s exp(.(V(x) - E)). (3.7) 

The correlation functions of the variables {V{x) — E)^, and therefore the moments of 
the effective Lyapunov exponent, can again be expressed as explicit multiple integrals. 
We will only need the expression of the mean effective Lyapunov exponent, which reads 

The rescalings z — x/icr, u = sE show that the above result scales as 
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Figure 4: Plot of the variance of the effective Lyapunov exponent for a Brownian potential 
{H = 1/2) with E = 0, against the logarithm of the scaling variable X = NA^. Full curve: 
analytical prediction (2.21). Dashed line: limiting value (A.27): varGi/2 = (tt^ — 4)/16 = 
0.366850. Symbols: numerical data for A = 0.1 and A = 0.2. 



in terms of the length ratio 

X = ^, (3.10) 

^cr 

where the energy-dependent crossover length 

4r = (^/A)^/^ (3.11) 

is the typical size of the first allowed zone containing the origin. 
3.2 Short samples 

For rather small samples, such that L <S i^r, i-e., X -C 1, the contour integral in (3.9) is 
dominated by a saddle point at Uc ~ z~'^^ » 1. We are thus left with the exponentially 
small estimate 

(7(L)) ~ exp (-^) • (3.12) 

This expression again scales as the probability that V{L) is equal to so that a forbidden 
zone has just been entered. 

For a weak but finite disorder strength A, the mean effective Lyapunov exponent 
again has a finite limit of order A^ for X <^ 1. The onset of superlocalization again takes 
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place for a system size such that 



X^|lnArV(2i.)«i, L^[^-^^^y\ (3.13) 

3.3 Long samples 

The superlocahzation phenomenon is fuUy developed in the converse regime of long sam- 
ples (L 3> 4r, i-e., X 3> 1). The contour integral in (3.9) is dominated by the square-root 
branch cut at s = 0. We thus obtain a power-law growth of the mean effective Lyapunov 
exponent: 

irrespective of the energy E. 

The full distribution of the effective Lyapunov exponent again simplifies in the limit 
of a long sample. Indeed, for x ^ £cr, E is most often negligible with respect to V{x). 
Equation (3.6) therefore simphfies to 

l{L)^y^'^dx^{V{x)U. (3.15) 

The rescaling of x by L yields the following asymptotic form of the effective Lyapunov 
exponent of long samples: 

7(L) ^ E^/^ Yh ^ E^/^ {L/£,,f/^ Yh ^ (AL^) Yh, (3.16) 

where the fluctuating factor Yh is another non-linear functional of the normalized frac- 
tional Brownian motion: 

Yh = f^dx^{v{x))+. (3.17) 

This functional is investigated in the Appendix. It has a universal distribution which only 
depends on the Hurst exponent H . Its flrst two moments are evaluated for all values of 
H [see respectively (A. 6) and (A. 22)], as well as the full distribution of Yi in the ballistic 
limit. Figure 5 illustrates the dependence of the dimensionless ratio 

W„ ^ m (3.18) 

on the exponent H . Figure 6 shows a plot of the probability density of Yi (ballistic case: 
exact expression (A. 33)) and of Y'1/2 (Brownian case), the latter being measured from an 
ensemble of numerically generated, and suitably rescaled, long random walks. The delta 
peak at F = 0, which is present in the ballistic limit {H = 1), gets smeared out into a 
continuous density which diverges as y — > for any H < 1. 
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Figure 5: Plot of the moment ratio Wh, defined in (3.18), of the quantity Yjj entering the 
result (3.16), evaluated from the exact expression (A. 23), against the Hurst exponent H. Sym- 
bols: results (A. 25), (A. 28), and (A. 30) for the particular cases H = 0, 1/2, and 1. 



The present work is devoted to the specific features of localization in a weak self-affine 
disordered potential in one dimension, whose fluctuations grow with distance with a 
positive Hurst exponent H. The most salient feature of such non-stationary potentials is 
that the effective disorder strength depends on the spatial scale. Even in the regime of 
most physical interest, where the fluctuations of the disordered potential are small at the 
microscopic scale (A <^ 1), the effective disorder becomes strong beyond a well-defined 
crossover length Ncr or £cr, which diverges as A~^^^ . Samples much longer than this 
crossover length typically consist of an alternation of classically allowed zones, where the 
growth rate of a generic wavefunction is small and proportional to the square disorder 
strength A^, and of forbidden zones, where this growth rate is large. Typical eigenstates 
on such samples turn out to be superlocalized [12, 13]. 

In this paper we present a parallel analytical investigation of the statistics of su- 
perlocalized eigenstates, for both the discrete tight-binding model and the continuum 
Schrodinger equation. First of all, the clear separation of the growth rates of wave- 
functions in allowed and in forbidden zones in the weak-disorder regime is shown to 
allow for a semi-classical treatment of the superlocalization phenomenon, by means of 
the W.K.B. formalism. The key quantity is the size-dependent effective Lyapunov expo- 
nent, which is closely related to the transmission across the sample, and therefore to its 
zero-temperature conductance by the Landauer formula. 



4 Discussion 
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Figure 6: Plot of the probability density of the random variable Yh entering the result (3.16). 
Thin full line: continuous component (delta peak at y = omitted) of the exact result (A. 33) 
for H = 1 (ballistic limit). Line with symbols (showing histogram bins): numerical result for 
H = 1/2 (Brownian case). 

For the discrete tight- binding model, investigated in Section 2, the effective Lyapunov 
exponent behaves as 

-fN-H ln(iV/iV„) +Gh^ ln(AA^-^) + Gh (4.1) 

in the fully supcrlocalized regime of long samples (A^ ^ Ncr) [see (2.26)]. The fluctuating 
part Off is an additive term of order unity. It is given by the functional (2.27), whose limit 
distribution only depends on the Hurst exponent H. These results translate as follows in 
terms of physical quantities. One of the most characteristic features of superlocalization 
is the fall-off of the effective locahzation length as a function of N. This decay is 
however only logarithmic, as ~ ^/[Hln[N/Nci)]. As far as the transmission Tjv is 
concerned (or equivalently the conductance qn), the relevant quantity is its logarithm, as 
in most localization problems. Its mean grows as (InT/v) ^ —2HNl'n.[N/Ncr), whereas its 
variance grows as varlnTV ~ 4 var G^A^^. The logarithm of the transmission is therefore 
marginally self-averaging, as its relative variance var lnTjv/(lnTjv)^ ~ l/[ln(A^/A^cr)]^ falls 
off logarithmically. 

For the continuum Schrodinger equation, investigated in Section 3, the effective Lya- 
punov exponent behaves as 

7(L) ^ E'/' {Lji^YI^ Yh ~ (AL^) (4.2) 

[see (3.16)] in the fully supcrlocalized regime. The fluctuating part Y^ of this result, given 
by the functional (3.17), is now involved as a multiplicative factor. This multiplicative law 
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generates more pronounced superlocalization effects, and especially stronger fluctuations 
in physical quantities. The effective localization length ^{L) ^ l/(y^A^/^L^/^) now 
decays as a power of the sample length, with a fluctuating amplitude. As far as the 
transmission is concerned, its mean grows as (In T(L)) ~ — A^/^L^"'"'^/^, i.e., more rapidly 
than linearly with the sample length. The logarithm of the transmission is now fully non- 
self-averaging, as its relative variance var In T(L)/ (In T(L))^ — > Wh — 1 has a non-trivial 
limit for L ^ £cr- 

The effective Lyapunov exponent and related physical quantities therefore exhibit 
two different kinds of statistics in the discrete and in the continuum formalism. The 
basic difference between both situations is already present in the case of a constant 
potential. It is indeed deeply rooted in the dispersion relations of the models. For the 
tight-binding model with a constant site potential V outside the band, the growth rate t of 
wavefunctions is such that E — V — 2 cosh t, and therefore only diverges logarithmically 
as t ~ ln|y| for large V. For the continuum Schrodinger equation with a constant 
potential V, the growth rate now obeys V — E = t'^, and diverges as t ~ ^/V. Replacing 
in the above dispersion estimates t by 7, and V by the product of a random number of 
order unity by Kms ~ AN^ or AL^ , leads to the correct qualitative forms of the scaling 
equations (4.1) and (4.2). Finally, this line of reasoning also shows that the results for the 
Schrodinger equation cannot be recovered from those of the tight-binding model by going 
to the continuum limit, because this limit corresponds to t — > 0, whereas superlocalization 
is essentially due to large values of t. Table 1 summarizes the discussion. 



Quantity 


tight-binding model 
N > A^cr sites 


Schrodinger equation 
length L » 4r 


Anderson 
localization 


-(InT) 
varlnr/(lnr)2 


l/hi{N/N,,) 
N\n{N/N^,) 
l/[ln(Ar/7V,,)]2 


(AL^)-V2 
constant 


1/A2 
A^L 



Table 1: Scaling behavior of various quantities in the regime of fully developed superlocal- 
ization: mean effective loealization length, mean and relative variance of the logarithm of the 
transmission (conductance). Only the qualitative scaling behavior is given (dependence on the 
sample size N 01 L and on the disorder strength A). The emphasis is put on the difference 
between the discrete tight-binding model (additive fluctuation in the effective Lyapunov expo- 
nent) and the continuum Schrodinger equation (multiplicative fluctuation). The well-known 
results for conventional Anderson localization are recalled for comparison. 

Let us close up with a word on the generality of our results. For any weak self-affine 
(not necessarily Gaussian) random potential obeying power laws of the form (1.2) or (3.2), 
with any positive Hurst exponent H, both the functional form of the results (4.1), (4.2) 
and the scaling laws recalled in Table 1 can be shown to still hold true. The fluctuating 
parts Gh and Yh arc, however, only given by (2.27) and (3.17) for the Gaussian self- 
affine potentials associated with fractional Brownian motion. Their distribution involve 
in general further details of the statistics of the potentials. 
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Appendix. The functionals Gh and Yh 

This Appendix is devoted to the random quantities Gh and Yh, which respectively enter 
our key results (2.26) and (3.16). They are defined as the following non-hnear functionals 

G^^ = ^'da:ln|^;(a:)|, Y„ ^ j\x ^ {v{x))+, (A.l) 

of the normalized fractional Brownian motion v{x) with Hurst exponent < if < 1 on 
< X < 1, so that v{Qi) = and 

{{v{x)-v{y)f) = \x-y\^^. (A.2) 

First moments 

The means (first moments) 

{Gh)^ f^dx{\n\v{x)\), {Yh) = J^' dx ^^(^;(x))+) , (A.3) 

can be readily evaluated. Indeed, v{x) is a Gaussian with variance (v(x)^) = x"^^ , so 
that its probability density reads 

Elementary integrals yield 

(ln|^(x)|) = Hlnx-^^, {^{v{x)u) = ^V^^^ ^ ft) ^""^^ ' ^^"^^ 
so that 



Second moments 



The second moments 

{G%) = 2 £ dx £ dy (In |t;(a;)| In \v{y)\), {Y^) = 2 J^' dx £ dy y{v{x)Uv{y))^'j , 

(A.7) 

now involve two-point observables, expressed in terms of the correlated Gaussian variables 
v{x) and v{y), such that 

{vixf) = x^^, {viy f) = y^^, {vix)v{y)) = i (x^^ + - {y - xf"") . (A.8) 
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It is convenient to use polar co-ordinates in the v{x), v{y) plane. For x < y, we set 

z = -<l, (A.9) 

y 

and 

r 

v{x)=rcos9, v{y) = — sm{9 + a), (A. 10) 

with 

{v(x)v(y)) 1 + ^2^- (1-^)2^ , , , , 

sing^ ^ )f ^— (0<a<n 2). A.ll 

^{vixf) (viyY) 

This parametrization ensures that the angle 9 is uniformly distributed between and 27r, 
whereas the radial variable r has a probability density 

^^W = ^exp(-^). (A.12) 

The calculation of (In \v{x)\ In li'(y)l) splits into the radial integrals 

(lnr) = ^ ^ , ((lnr)2)e = -, (A.13) 

where (• • •)c denotes the connected part, and the angular integrals 

(ln|cos^|) = (ln|sin^|) = -ln2, (ln|cos(^)| In |sin(^ + a)|)c = ^-7^- (A.14) 
The derivation of the second result requires the use of the Fourier series 

n=l n=l -"-^ 

where the latter equality holds for |^| < tt. We thus obtain the following simple expression 

(ln|^;(a;)|ln|^;(t/)|)e=y (A.16) 
for the connected correlation function, which leads us to the result 

■Gh^{GI),^\ [\^dz, (A.17) 
Jo 

where a is defined in terms of z by (A.ll). 

Similarly, the calculation of (^^ {v{x))^{v{y))^'^ splits into the radial integral 



van 



(r) = ^|x^ (A.18) 
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and the angular integral 

(A.19) 

(ot 7r\ 
e + - - -j = A;sin0, (A.20) 



where we have set 
with 



and where E(A;) and K(A;) are complete elliptic integrals. Finally, 

^^^^ ^ ¥T2 /I ll (^^^^ ~ ~ ^'^^^^^) ^^-^^^ 

so that 

where /c is defined in terms of z by (A.21). 

The general expressions (A. 17) and (A. 23) simplify in the following particular cases. 



Ultraslow limit {H = 0) 

In this limit we have 

irrespective of hence 



varGo = — = 0.137078, 
72 



= — , (A.24) 
6 2' ^ ^ 



(A.25) 



Brownian case {H — 1/2) 

In this case we have 

sin q; = 2A;2 - 1 = V^, (A.26) 
so that a direct integration yields 

-4 

varGi/a^ = 0.366850. (A.27) 

In order to evaluate Wi/2-i it is more convenient to go back to the middle expression 
of (A.19), and to integrate first over k, then over (j). We thus obtain after some algebra 

= r(|)' (50 - 3^2 In (l + V2)) = 1.895949. (A.28) 
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Ballistic limit (H ^ 1) 

In this case we have 

a = |, A; = 1, (A.29) 

irrespective of z, hence 

varCi = — = 1.233701, Wi = -^t(\Y = 2.360681. (A.30) 

The constancy of a = 7r/2 means that v{x) and v{y) are fully correlated. Their correlation 
functions (A. 8) indeed saturate the Cauchy-Schwarz inequality. 

The whole fractional Brownian process degenerates into a one- variable problem in the 
ballistic limit. We have indeed 

v{x) = xw, (A. 31) 

where w = t'(l) is a Gaussian variable such that (w^) = 1, hence 

Gi = In \w\ - 1, Yi = '^y/w^. (A.32) 

By performing the appropriate changes of variables on the Gaussian law of we get the 
following explicit expressions for the probability densities: 

P(GO = exp (G, + 1 - le^(^^+^)) , P{Y,) = lS{Y,) + ^ exp . 

(A.33) 
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